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Summary. We conduct 2D axisymmetrical hydrodynamical simulations to inves- 
tigate the interaction of a collimated fast wind (CFW; wide jets) with a spherical 
AGB wind. The code includes radiative cooling. We find that the shape of the plan- 
etary nebula (PN) is sensitive to the exact mass loss history of the AGB wind, and 
the opening angle of the CFW. Some typical PN morphologies are obtained, but 
many other observed morphologies seem to require more ingredients than what we 
assume in our present simulations, e.g., equatorial AGB wind, and ionization and 
fast wind during the PN phase. The hot bipolar bubble formed by the jets is an 
X-ray source. 
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1 Numerical Simulations 

Our simulations were performed using Virginia Hydrodynamics-I (VH-1), a 
high resohition multidimensional astrophysical hydrodynamics code devel- 
oped by John Blondin and co-workers (Blondin et al. 1990; Stevens et al., 
1992; Blondin 1994). We have added radiative cooling to the code at all 
temperatures T > T„iin, where the gas temperature is bound from below 
at Tniin — 300 — 1000 K. Radiative cooling is carefully treated near contact 
discontinuities, where large temperature gradients can lead to unphysical re- 
sults. 

We simulate axisymmetrical morphologies. This allows us to use axisym- 
metrical grid, and to simulate one quarter of the meridional plane. There are 
208 grid points in the azimuthal {9) direction of this one quarter and 208 grid 
points in the radial direction. The radial size of the grid points increases with 
radius. In these simulations the grid extends from 10^^ cm to 4 x 10^^ cm. 

Before the jet is launched [t — Q), the grid is filled with slow wind having 
a speed vi = 10 km s~^ and mass loss rate Mi ~ 1O~^M0 yr^^, with small 
variations between the runs. We launch a collimated fast wind from the first 
20 zones attached to the inner boundary of the grid. The jet is uniformly 
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ejected within an angle (half opening angle) a (0 < 9 < a). For numerical 
reasons a weak slow wind is injected in the sector a < 9 < 90° {9 — along 
the symmetry axis - vertical in the figures here). 

2 Results and Summary 

We try to explain PN shaping and formation processes by simulating jets 
interacting with the AGB wind. The parameter space for these types of flows 
is huge. We run more than 200 simulations, in which we tried many values of 
the CFW (jets) mass loss rate, velocity, and opening angle. 

Another parameter that was used is the mass loss rate history of the AGB 
wind. For example, in a model presented here (figure [2]), we assume that some 
period At before the beginning of the jet-launching phase the slow wind mass 
loss rate increased by some factor fc; e.g., in the run presented in figure [2] 
At = 950 yr and fc = 6. We found that the mass loss history of the AGB wind 
plays a significant role in the shaping process. In some runs the AGB wind 
initial temperature was set to 10, 000 K instead of 1, 000 K, with a noticeable 
influences on the final shape of the PN (compare the two panels in figure 
[2]). Taking the AGB wind to be at 10,000 K is relevant to system where 
the accreting companion is a WD. If the accretion rate is high enough the 
WD sustains a constant nuclear burning, and it is very hot and luminous, as 
super-soft X-ray sources (e.g., Starrfield et al. 2004). Such an accreting WD 
might maintain the entire AGB wind ionized. 

We were able to reproduce some physical properties of PNs, but not all 
desirable properties. It is very difficult, and might be impossible, to get the 
entire range of properties with our limited numerical code. We must add 
more ingredients to the code. Examples are pulsed jets, precessing jets, dense 
equatorial outfiow and ionization fronts at later times. 

We presented results of four numerical runs that migh match real PNs. 
Our main results are as follows: 

1. We show that different PNs morphologies could be strongly dependent on 
the AGB wind mass loss rate history. 

2. We have an evident that the final morphology is very sensitive also to the 
assumed initial temperature. For Ti — 1000 K and Ti = 10000 K we get 
two different evolutions of the planetary nebula (figure [2|). 

3. The dense finger along the symmetry axis (in figure [21 top) result from 
instabilities. In reality, it would be more extended and not so narrow. It 
is forced to the symmetry axis by the numerical code. 

4. In one case we obtain low density finger (in figure [2]) protruding to the 
upper right form a torus. However, it is also a result of an instability, 
and we expect that in real systems there will be several such fingers. In 
particular, if there is a departure from axisymmetry due to the orbital 
motion, we would expect the fingers to be similar but not identical in the 
two sides of the equatorial plane. 
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Fig. 1. Top: Density plot ai t — 2000 yr (after jet launching starts) for a case 
where the slow wind mass loss rate has increased by a factor of five to A4i = 
5 X 1O~^M0 yr~^ 950 yr before jet starts. The vertical axis is the symmetry axis 
and the horizontal axis is in the equatorial plane. Slow spherical wind: at t = 
vi = 10kms-\ Ml = 5x10-^ Mq yr"^ for r < 3x10^^ cm, and Mi = IO'^Mq yr"^ 
for r > 3 X 10^^ cm. The jet has a half opening angle a — 20°, V2 — 600 km s~^, 
and M2 = 2 X lO~**Af0 yr~^. Arrows indicate flow direction: v > 200 km s~^ (long 
arrow), 20 < u < 200 km s~^ (medium arrow), and v < 20 km s~^ (short arrow). 
This general case can account for PNs such as NGC 6886 or NGC 3918. Bottom: 
Density plot at t — 530 yr for another case. Slow spherical wind: vi = 10 km s~^. 
Ml = 5 X 1O"®M0 yr"^ The jet has a half opening angle a = 40° , U2 = 600 km s"\ 
and M2 = 2.5 x 10"'^ Mq yr"^ This general case can account for some morphological 
features in He 2-104, such as the instabilities in the outer lobes. 
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Fig. 2. Top: Density plot ai t = 1300 yr for a case with a slow wind of vi — 
10 km s^^, and mass loss rate of Mi = 3 x lO^^A/0 yr~^ for r < 3 x 10^'' cm, and 
Ml = 5 X 10"'' Mq yr"^ for r > 3 x 10^*^ cm at t = 0. The jet has a half opening 
angle a — 30°, 112 ~ 600 km s~^, and A42 = 1O~^M0 yr~^. The initial temperature 
was 1000 K, and the temperature is limited from below at 1000 K. Bottom: Density 
plot at t = 1560 yr for a case with a slow wind of vi = 10 km s~^, and mass loss 
rate of Mi = 3 x lO-'^M© yr~^ for r < 3 x lO" cm, and Mi = 5 x 10~*A'/q yr"^ for 
r > 3 X 10^'' cm at t = 0. The jet has a half opening angle a = 30°, V2 ~ 600 km s~^ , 
and M2 = 1O~^M0 yr"'^. The initial temperature was 10000 K, and the temperature 
is limited from below at 1000 K. The general shape is similar to MZ 3. The arrows 
are just as in figure [T] Note that top and bottom are the same beside the initial 
temperature. 
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5. In all cases that we run we saw that there is a strong dependence on the 
half opening angle of the jet. This will be the focus of a future paper. 
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